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Abstract 

We discuss some experimental signatures associated with the topological structures of unconven- 
tional superconductor order parameters of form drc2_y2+ix, where x = s,px±Py, or dxy Specifically, 
we study the topological surface states on the (110) and equivalent surfaces of such superconduc- 
tors which are observable in Andreev tunneling experiments, as well as evaluate the magnetic flux 
trapped in superconducting rings of such superconductors with multiple grain-boundary Josephson 
junctions. Previous experiments are examined and several new experiments suggested. 
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Definitive experiments relying on topological differences between different superconduc- 
tors to reveal their symmetry are of interest in themselves. The motivation for such investiga- 
tions in the Cuprates is especially strong. It has been suggested that the pseudogap regime 
in the phase diagram of the Cuprates ends at a Quantum critical point (QCP) and that 
the pseudogap regime represents an elusive time- reversal violating symmetry^. Symmetry- 
based proposals^i for experiments in the normal state using circularly polarized angle resolved 
photoemission experiments have been carried out. They show results^*^ consistent with a 
Time-reversal violating phase in the normal and superconducting phases below a critical 
density. If this is true, the superconducting order parameter for particle density below the 
QCP should also be a Time-reversal violating. 

The two time-reversal violating states suggested for the underdoped Cuprates have a 
current pattern sketched in fig. (1) of Reference(^. State I (magnetic group in tetragonal 
crystals 4/mmm) preserves symmetry on reflections in planes normal to the x — y plane 
going along the x and the y -axes but violates it on reflections in the x±?/ planes. It has two 
domains, specified by the sign of the phase or equivalently by the direction of the current 
in any of the four plaquettes in a unit-cell. State II (magnetic group mmm) has current 
pattern which violates reflection symmetry in planes going through the x and the y -axes. It 
has four domains two due to time-reversal and two because reflection symmetry is violated 
either over planes going through the {x + y) or the {x — y) -axes. This state breaks Inversion 
symmetry as well while State I does not. This is the state consistent with the experiments 
of Kaminski et al.— . Considerations of these symmetries^i^ and of domains is important for 
the design of the experiments. 

Suppose such states were to undergo continuous transitions to superconductors of dom- 
inant dx2-y2 symmetry at Tc. Then a time-reversal violating component would be linearly 
admixed which by continuity preserves the reflection symmetry of the Time-reversal violat- 
ing normal state. For State I, such an admixed state is uniquely of the d^y variety giving 
a dx2_y2 ± id^y superconductor; The ± sign signifies the two time- reversed domains. For 
state II, the admixed states possible are either of the Px+y oi Px-y variety, depending on the 
domain, giving a dx2-y2 ± i{px ^Py) superconductor. (Due to lack of Inversion symmetry a 
triplet superconducting state is linearly admixed to the singlet). Kaur and AgterbergS have 
already come to these conclusions on formal group-theoretical grounds. 

Suggestions for time-reversal breaking through admixture of the is or idxy type near sur- 
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faces due to de-pairing of the principal d^-^-y-i order parameter and their effect on Andreev 
reflection experiments have already been studied^. Our principal interest here is not such 
surface induced time-reversal breaking which would set in at a lower temperature than the 
bulk Tc but time-reversal breaking in the bulk, possibly existing only in cuprate supercon- 
ductors for hole or electron doping less than at the quantum critical point. There have 
been proposals for time-reversal violation in the superconducting state alone from differ- 
ent considerations^^ which all preserve center of inversion. We consider below three types 
of Time-reversal violating superconducting states: d^-^-yi + ix] [x = s, d^y or Px±y]- The 
imaginary component will be considered to be small compared to the real part. Where the 
context is clear, let us call them d + is,d + id,d + ip. For the first two cases the symmetry 
of the Andreev effect due to bulk time-reversal breaking is the same as due to surface in- 
duced time-reversal breaking already studied^ but the magnitude of the effect can be quite 
different. 

To determine the precise form of order parameter experiments that can detect the pres- 
ence of the small additional component unambiguously must be designed. We discuss here 
two kinds of plausible experiments. The experimental signatures studied are consequences 
of the topological structures of the order parameters and are robust to small perturbations. 
The first experiment discussed is a tunnelling experiment on the (110) (and equivalent) 
surfaces of the high- Tc cuprates of the Andreev type. For the x = s, dxy, the symmetry 
of the Andreev effect due to bulk time-reversal breaking is the same as due to surface in- 
duced time-reversal breaking already studied^ but the magnitude of the effect can be quite 
different. 

A. surface states on (110) surfaces 

The presence of surface state with energy e = on the (110) (and equivalent) surfaces 
of (pure) (i^2_j^2-wave superconductors was first pointed out by Hu^°. The existence of this 
zero energy surface state is a direct consequence of the topology of the superconducting 
order parameter. We extend Hu's analysis here to study surface states in superconductors 
with a small deviation from the pure d^^^-y^ symmetry. Following Hu's notation, we rotate 
our co-ordinate axis by 45° on the x — y plane such that the x— and y— axes in the new 
coordinate system are the (110) and (110) crystal axes of the cuprates, respectively. In 
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this notation, the order parameter of a pure (ia.2_y2-wave superconductor has the property 
A{kx, ky, kz) = —A{—kx,ky,kzY-. The electronic excitations of an inhomogeneous super- 
conductor are determined by the Bogohubov-de Gennes equations, 

€nUn{x) = HoUn{x) + j c/x A (f , f )f „ (f ) (l) 
^nVn{x) = -HoVn{x) + j (ifA(f , f 

where Ho = ~2m^i — + v{x), v{x) is the electron crystal potential energy and A{x,x') 
is a general, non-local superconductor order parameter. A superconductor at a; > with 
surface at x = with thickness d can be modelled by a superconductor order parameter of 
form 

A{k,x) ~ Ao{k)e{x), 

with boundary condition that the electronic wavefunctions all vanish at x = —d^^. Ao{k) is 
the wave-vector-dependent order-parameter of the bulk superconductor. The Bogoliubov- 
de Gennes equations for electronic states close to the Fermi surface can be solved in the 
WKBJ approximationi«^ . In particular, surface states with u{v)(x) ~ Q^ikpyy+kF^z)^-"/x ^^^^ 
determined by the following eigenvalue equation derived by Hu^-, 

W ) [, - iy^\K4^ J ^' 

where |e| < |A+(_)| is the surface state energy, A+ = Ao^kpx, kpy, kpz) and A_ = 
Ao( —kpx, kpy, kpz). kp — {kpx, kpy^ kpz) is a wave- vector on the Fermi surface. 

Notice that e depends on kpx,kpy implicitly through the momentum dependence of the 
gap function AQ{k). For pure d-wave superconductors, A^ = — A_ and the eigenvalue 
equation has a particular solution e = that is independent of momentum kp and surface 
thickness d. This is the zero energy mid-gap state discovered by Hu-*^- . Note that other 
surface states may also exist. However they always exist in pairs with energies ±e because 
of particle-hole symmetry in the eigenvalue equation. The mid-gap states can be probed by 
tunnelling experiment across the (110) surface which measures the density of states directly. 
A "zero-bias conductance peak" is expected to be observed with a particle-hole symmetric 
tunnelling spectrum. 

Now consider superconductors with a small non-d-wave order parameter component. In 
this case the relation A_|_ = — A_ is modified. To begin with, let A be real with A_|_ = A + s 
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and A_ = —A + s, where s is independent of k, corresponding to the allowed symmetry in 
an orthorhombic crystal with Time-reversal preserved. The case of complex gap function 
which breaks time-reversal symmetry will soon be discussed. By direct inspection of Eq. 
(j21) we see that the e = surface states that exist for all momenta kpy, kpz for pure d-wave 
superconductors, survive in the presence of the small s-component, as long as /S.{kpy, kp^) > 
s. Other surface states may also exist in pairs with energies ±e because of particle-hole 
symmetry. Since the surface states energies must satisfy |e| < A{kpy, kpx) — s, the main 
effect of mixing a small, real non-d-wave component is to reduce the number of available 
surface states that are allowed in pure d-wave superconductors. Consequently, tunnelling 
experiment cannot clearly distinguish pure d-wave superconductors from superconductors 
with a small, real, non-d-wave component. 

The situation changes if we consider superconductors with complex order parameter 
satisfying A+ = — Al. In this case the eigenvalue equation becomes 



different domains of time-reversal or direction of flow of internal currents. The particle-hole 
symmetry that exists in pure d-wave superconductors is lost because of broken time-reversal 
symmetry and e = state is no longer a solution to the equation. Eq. (3) has the formal 
solution 



Writing e = |Ao|sm6' it is easy to see that 6 = (p in the limit d = and the surface state 
energy measures the imaginary part of the gap function |Ao|sm0 directly, e can be either 
positive or negative, depending on the sign of (p. For d ^ the surface states do not measure 
the imaginary part of gap function directly and both positive and negative energy surface 
states occur in general. 

Let us apply the above result to see if x with different symmetries can be distinguished in 
dx2-y2 +ix superconductors. In figure (1) we show the energy gap structure on the kj> — ky 
plane for superconductors with such order parameters, for x = s,p,dxy First we consider 
dx2-y2 + is and d^^-y^ + idxy superconductors. Both superconductors have the property that 
A —A* upon refiection on the (110) and (llO) planes, implying absence of zero energy 





is defined through A_|_(_) = |Ao|e^*^ The two signs in front of correspond to the two 
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FIG. 1: order parameter structure on the kx — ky plane for superconductors with order parameter 
dx2_y2 + i5 for 6 = s,p,dxy Notice that there is another configuration for the d + ip state where 
left and right region interchanges 

surface states and appearance of particle-hole asymmetry in the surface state spectrum for 
any fixed kpy However, particle- hole asymmetry vanishes in the tunneling spectrum if the 
tunneling experiment samples all kpy vectors equally, since A —A* when ky —>■ —ky, and 
the broken time-reversal symmetries in lower and upper ky planes are in exactly opposite 
domains. After adding up the contributions from all kpy wave vectors what should be 
observed is a particle-hole symmetric spectrum with the splitting of the zero-energy bias 
peak into two peaks at finite energy, as pointed out in a number of earlier worka^^i^. 

The situation is however, different for d^^-yi + i{px ± py) superconductors. While the 
reflection symmetry of order parameter upon the (110) and (110) planes are the same for 
the d + is and d + id states it is different for the d^^-y^ + i{px ± Py) state. A —>■ —A* upon 
reflection on the (110) plane but A —A on the (110) plane or vice versa. Furthermore, the 
domain of broken time-reversal symmetry remains the same (A —A) when ky — > —ky. 
As a result, the e = surface states will be absent on the (110) (or (110)) surface of 
d + ip superconductors and the corresponding tunnelling spectrum, after sampling over all 
kpy vectors, will remains particle-hole asymmetric, whereas e = surface states exist on 
the (110) (or (110)) surface with particle-hole symmetric tunnelling spectrum as in pure 
d-wave superconductors. The existence of different tunnelling spectrums on the (110) and 
(110) surfaces with one of them being particle-hole asymmetric is a very strong indication 
of d + ip superconductors. But experiments must be done in monodomain samples where 
(110) and (110) can be distinguished. There is the further problem that even in a single 
crystal, where for a given sign of time- reversal, the two different possible reflection domains 
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{px ± Py) may occur. For orthorhombic symmetry, this may be less likely. There is also the 
possibility that the boundary or tunnelling current may pin the orientation of the p order 
parameter in experiment and only one of two possible surface behaviors is observable. Thus 
if a strong asymmetry is observed in conductance with a peak predominantly for > or 
V < 0, a. bulk d + ip state may be suspected. 

For X = is, dxy, one must also consider the possible surface induced pair-breaking^. This 
can be distinguished by the temperature where it sets in, the fact that there is no obvious 
reason for a composition dependence for it and that its amplitude is expected to be smaller 
than the bulk effect if the thickness of the surface layer is small O (thickness of surface 
layer /length scale of bound state). Otherwise experiments which are sensitive to the bulk 
(but not surface) state order parameter have to be considered. We shall examine such 
experiments in section C. 

B. Existing Andreev reflection Experiments 

Given the above information let us look at the existing tunneling experiments to see 
what conclusions can be reached. Early Andreev reflections experiments on samples of 123- 
cuprate^^ with tunneling in the (110) and (110) directions have revealed symmetric peaks 
at zero voltage which split on applying a magnetic field as predicted^. These confirmed that 
the pairing is of d^i-yi symmetry. Subsequently two groups^^^ have reported experiments 
which in zero magnetic field show two peaks at ±SV much less than the superconducting 
gap. The two peaks have in general different heights. The authors correctly conclude 
that a peak not at zero-voltage is an indication of time-reversal violating superconducting 
order parameter. Particularly interesting is the fact that in a given tunneling junction, 
the asymmetry in the amplitude of the peaks increases in annealing. Asymmetry of about 
a factor of 4 has been reported^^. Unambiguous further experimental confirmation of the 
asymmetry on monodomain samples are therefore important because one of the exact results 
of our analysis is that the absence of particle-hole symmetry in the peaks in the tunneling 
spectrum indicates a bulk time-reversal violating d + ip state. Moreover, differences of 
opinioriF^ exist about the surface composition of the samples in which the bound state(s) 
exists at zero-voltage or at finite voltage. At this point one cannot be sure where the region 
of the peak at nonzero V lies in the phase diagram of 123-cuprates. 
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FIG. 2: arrangement of the two superconducting grains in the two-grain ring. Notice that the 
(100) axis of the d + i6-wave superconductor is 22.5" off the x-direction. 

C. Superconducting rings with multiple grain boundary Josephson junctions 

Next consider flux-quantization experiment on superconducting rings with multiple grain 
boundary Josephson junction o^^i^^i^^ . With suitable arrangement of the position of the grain 
boundaries these experiments can in principle be used to distinguish between different broken 
time-reveral symmetry states. Start with a simple ring made up of two mono-crystal grains 
as in the Wollman et al. experiments^^. One grain is a simple s-wave superconductor, and 
the other is the d + ix-wave superconductor of interests. The arrangement of the two grains 
is shown in figure (3). 

The (100) axis of the d + ix-wave superconductor is chosen to be 22.5° off one of the 
directions of the cut. Consider the flux quantization condition of the ring as a function 
of the cutting angle 6. From figure one we see that for 6 < 22.5° the phase of the d + ix 
superconductor remains constant over the whole region for all three choices of x. There is no 
mismatch in phase across the grain boundary, and consequently the flux trapped in the ring 
is integral multiple of magnetic flux quantum $o = h/2e, as in usual s-wave superconducting 
rings. For 22.5° < 9 < 67.5°, the phase of the dx2_y2 + idxy superconductor differs by 20 on 
the two surfaces but remains the same for the other two cases. Consequently the flux trapped 
in the ring is {n±^)^o for the d^^^yi+idxy superconductor, and remains equal to n$o for the 
other two kinds of superconductors. A different regime is reached for 67.5° < 9 < 112.5°. In 
this regime, the phase difference between the two surfaces of the d^^-yi +idxy superconductor 
is precisely vr, whereas it is 7r±20 for dx'2-y2+is superconductor, and can be either tt or 7r±20 
for the dr^2_y2 + ip superconductor, depending on the orientation of the p-order parameter. 
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FIG. 3: arrangement of the three superconducting grains in the experiment by Tsuei et.al. 

Provided that the phase is not too smaU, the different behaviors of magnetic flux trapped 
in these two regions of ^- value can distinguish d + id superconductor from the other two 
types. The presence of two different possible values of trapped flux, if observable, can also 
distinguish d + ip superconductor from the rest. However the magnetic flux may pin the 
orientation of the p-order parameter, or if only one of the two {px ^Py) arc present, only one 
type of behavior is observable. The d + ip wave superconductors can be distinguish from 
the rest if we further increase 6 to 157.5° < 6 < 202.5°. In this regime, the phases on the 
two surfaces of the d + ip superconductor differ by 20 but remain the same for the other 
two cases. The trapping of non-integral flux quantum in this region is a clear signature for 
d + ip order parameter. 

Another interesting configuration that can distinguish d + is and d + id superconductors 
from the rest is to study Josephson effect across two 180° grains in which both grains are 
of the cuprate superconductors with surfaces pinned at the ±(110) directions. This is the 
direction which is gapless for pure d^2_y2-wave superconductor but are gapped for d + is 
and d + id superconductors. The existence of Josephson effect in this orientation provides 
unambiguous proof that the superconductor has a mixed order parameter. In this geometry, 
the d + ip wave superconductor presents a very interesting possibility in principle. Since 
d+ip wave superconductor is gapless along one of the (110) or (110) directions and gapped in 
the other, Josephson effect with phase shift near 7r/2±0 is to be expected in one orientation 
and no Josephson effect in the other. There are probably experimental limitations for such 
observations which we are not experts on, as well as the possibility that the orientation of 
the p order parameter may be pinned by the boundary. 
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Finally, let us study the situation of a 3-junction ring in the celebrated geometry first 
fabricated by Tsuei et.al}'^ . The ring is made by the superconductor of interests only, but 
with different crystal orientations in the three grains. On the — ky plane the three grains 
cover regions specified by angles 180° — 360°, 303.5° — 360°, and 26.5° — 150°, (0° is the 
(110) axis) respectively (see figure three). For pure c/-wave superconductor there is a phase 
difference of vr on the two surfaces of the second grain but the phases are the same on the 
two surfaces of other two grains. Consequently the flux trapped in the ring is precisely 
(n + i)$o. 

Repeating the analysis for the three different kinds of time-reversal violating supercon- 
ductors of interest here, we find that the flux trapped in the corresponding rings of d + is 
superconductor is given by F, = (ra + 1 ± ^)<l>o, whereas the flux trap in the d + ip super- 
conductor can be either Fp = (n + |)$o or (n + ^ ± ^)<l>o, depending on the orientation of 
the p-order parameter. The flux trap in the d + id superconductor is the same as d + is 
superconductor, Fd = (n + | ± ^)$o- However, there is some ambiguity because of the 
node in idxy order parameter along x and y axes. Note that with different cutting angles 
the value of the trapped flux will change. For example, if the three grains cover regions 
specified by angles 180° — 360°, 280° — 350°, and 60° — 170°, then the flux trapped in a pure 
d-wave superconductor ring will be n$o, whereas the d + ip superconductor will trap flux 
Fp = {n i: ^)$o; independent of the orientation of the p-order parameter. 

Recently, Tsuei et. al. have reported results of experiments on superconducting rings with 
essentially their original 3-junction ring geometry for three different cuprate systems^-^ for 
various dopings. They observed that the trapped flux is always equal to (n + |)<l>o within 
experimental accuracy of a few percent. Based on this observation, they conclude that the 
mixing of dx2-.y2 + is and dri.2-.y2 + idxy states, if exists, must be small in a large range of 
dopping in high-Tc cuprates. While we agree with their conclusion on the dx2-y2 + is and 
our analysis shows that the dx2-y2 + i{px ± Py) state is not ruled out 
by their results. We therefore suggest experiments with different cutting angles as discussed 
above. These may also be more precise since the dominant flux due to the dominant d-wave 
symmetry is nulled out in such experiments and only the corrections are visible. 
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